Perturbations of generic Kasner spacetimes and their stability 
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This article investigates the stability of a generic Kasner spacetime to linear perturbations, both 
at late and early times. It demonstrates that the perturbation of the Weyl tensor diverges at 
late time in all cases but in the particular one in which the Kasner spacetime is the product of 
a two-dimensional Milne spacetime and a two-dimensional Euclidean space. At early times, the 
perturbation of the Weyl tensor also diverges unless one imposes a condition on the perturbations 
so as to avoid the most divergent modes to be excited. 

PACS numbers: 98.80.-k,98.80.Cq 



I. INTRODUCTION 



The formalism of cosmological perturbations about a 
Bianchi I universe with a scalar field was only developed 
recently in Refs. [l|, and in Ref. 0] for the subcase of 
axisymmetric spacetimes. In such anisotropic inflation- 
ary models, the shear is always dominating at early time 
so that the universe behaves as a Kasner spacetime. It 
was realized that this anisotropic early era has an impor- 
tant signature since gravity waves are amplified during 
this era 0,13 • Indeed, this stage is usually short and this 
instability is transient. 

However this has led to question the stability of a pure 
Kasner universe [1]. The Kasner spacetimes [H are vac- 
uum solutions of the Einstein field equations. They de- 
scribe universes which are spatially homogeneous and Eu- 
clidean but with an anisotropic expansion. They play an 
important role in cosmology since they are a key struc- 
ture in the discussion of the dynamics of spatially ho- 
mogeneous spacetimes close to the singularity. Belinsky, 
Khalatnikov and Lifshitz 0-0 and Misner Q investi- 
gated the nature of the cosmological singularity by means 
of a Bianchi IX model, whose temporal behaviour toward 
the singularity was shown to be described by a sequence 
of anisotropic Kasner era. This has led to the mixmaster 
picture @ and the idea of the cosmic billiard that 
can offer a description of the geometry of the universe 
prior to infiation. 

The stability of this picture was already adressed in 
Ref. @ and later numerically in Refs. [TT| - [l3| . Only re- 
cently was it revisited Q in light of the recent develop- 
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ments concerning Bianchi I universe perturbations [ll-Q 
in the particular case of a Kasner universe with an axial 
symmetry, but no general study of the stability to linear 
perturbations of a Kasner spacetime with arbitrary ex- 
ponents has been performed yet. This is the goal of this 
article. 

We limit our analysis to the study of the stability of 
the Kasner spacetime to linear perturbations. Indeed, 
this can exhibit some instabilities but does not allow to 
demonstrate the general stability when no instability is 
exhibited at linear order in the perturbations. For in- 
stance, the stability of the Minkowski spacetime to lin- 
ear perturbation is an easy exercise while the general 
proof of stability is extremely challenging |14| . Also some 
spacetimes may be unstable to a particular class of per- 
turbations while stable in other cases. For instance, the 
Einstein static universe was shown to be unstable with 
respect to spatially homogeneous and isotropic pertur- 
bations [15[ while it was then understood that the issue 
was not that simple when Harrison [l^ demonstrated 
that all physical inhomogcncous modes arc oscillatory 
if the matter content was a radiation fiuid, while Gib- 
bons demonstrated the stability against conformal 
metric perturbations when the matter content was a fluid 
with sound speed larger that i. Finally [13, it was 
shown that the Einstein static spacetime is neutrally sta- 
ble against small inhomogcncous vector and tensor per- 
turbations and neutrally stable against adiabatic scalar 
density inhomogeneities as long as the condition exhib- 
ited by Gibbons holds, and unstable otherwise. Our anal- 
ysis is somehow simplifled by the fact that the Kasner 
spacetimes are empty so that we need no assumption on 
the matter content or on the type of perturbations con- 
sidered. Once the evolution of the linear perturbation 
is determined, in a gauge invariant formalism to ensure 
that no gauge mode can spoil the conclusion, we compute 
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the invariants of the metric and compare them to those 
of the background. Again, since the Kasner spacetime is 
empty, only the square of the Weyl tensor has to be con- 
sidered. This ensures that the conclusion is not affected 
by the choice of gauge or coordinates system. 

In this article, we consider the perturbation theory 
around generic Kasner spacetimes, as described in § 
in order to discuss their stability. First, in § IIIIl starting 
from the general formalism we developed in Rcfs. [l|, [^j, 
we show that for a vacuum spacetime only two pertur- 
bations are propagating, as expected. We then exhibit 
some asymptotic analytical solutions for the behaviour 
of the perturbations in § IIVI which allow to discuss the 
stability at late and early times in the generic case. In 
§ El we consider the axially symmetric case considered 
in Ref. Q , mostly in order to discuss the agreement with 
this previous analysis. 
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FIG. 1: Evolution of q; with vj (qi in solid line, §2 in dashed 
line and qa in dotted line). We can choose £ [-1, |-] so that 
Pi < P2 < P3. 



II. BACKGROUND SPACETIME 

A. Definitions 

We consider a Kasner spacetime, that is a Bianchi I 
universe with metric usually written in terms of the cos- 
mic time t as 



(2.1) 



for the vacuum, where S is the volume averaged scale 
factor and 'Jijit) the metric on the constant time hy- 
persurfaces. The vacuum Einstein equations imply that 
H = —3H^, from which we deduce that the averaged 
Hubble parameter is 



H : 



s 



1 



(2.2) 



It is easily integrated to deduce that S(t) = So{t/to)3 
and we can always set 5*0 = 1. 

The Kasner metric can then be shown to be of the form 



d.s^ 



-dt' 



E 



dx' 



with the coefHcients Pi satisfying the constraints 



that are fulfilled if 
2 . 



Pi 



■ sm Wi 



277. 



(2.3) 



(2.4) 



(2.5) 



zu varies in an interval of length 27r/3 and we always 
have pi < {p2,P3) if w G [—n/Q, 7r/2] but it is enough to 
let it vary in w G]7r/6, 7r/2[, which always ensures that 
Pi < ^2 < Ps (sec Fig. [IJ. The two limiting cases ru = ^ 



or w ^ ^ for which p.; = { — i, |, |} and pi ~ {0, 0, 1} re- 
spectively, correspond to space-times with an extra axial 
symmetry so that p2 = P3 or p2 = Pi- 

It is thus clear that the spatial metric is explicitely 
given by 



2<?, 



with q^=pi-^ (2.6) 



= diag 
so that 

i i 

The Qi are explicitely given by 



= - sniro;. 
We also define the useful parameterisation 



Qi = 2* ^ sinwi 



(2.7) 



(2.8) 



(2.9) 



such that Qi e [—1,1]. We define the conformal time ?/ 
from dt = Sdij, and the metric (j2.1[) can then be recast 
as 



ds' 



S^irj) [-drj''+j,,(rj)dx'dx^] 



(2.10) 



From the previous analysis, we have rj = |io(V^o)^ so 
that we set 



and also 



Siv) = with rjo ^to 



(2.11) 



(2.12) 
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We easily deduce that 

Mi 



7„- = diag 



7'-' ~ diag 



(2.13) 

We define the comoving Hubble parameter by H = S' / 
where a prime refers to a derivative with respect to the 
conformal time so that 



n = SH 



2^ 



The shear tensor is defined as cr, 
(7*Vy- = g\ = 0) so that 



(2.14) 



7- and is traceless 



^ = 3diag 



3qi 



and thus S*- = 3diag[(7i] so that 



6, 



(2.15) 



(2.16) 



and the Einstein equations imply that the shear satisfies 



2-Ha! 



0. 



\ 3/ ' '"J 
We also define the shear in cosmic time by 



ld7 



2 dt 



S 



(2.17) 



(2.18) 



B. Weyl Tensor 

The Kasner spacetime is a vacuum solution of Einstein 
equations and this implies that 



Ruv — 0, i? — 0, R 



■fii/pa 



(2.19) 



It is thus only characterized by its Weyl tensor, both at 
the background and perturbation level. The Weyl tensor 
has just two types of non-vanishing components which 
are the C'^'q^- and the C'°*^j.. The former are related to 

the components of the type C^\i (thanks to the traceless 
conditions) and to the electric part of the Weyl tensor. 
The latter are related to its magnetic part, and it vanishes 
for the background Kasner spacetime. The only non- 
vanishing components of the background Weyl tensor are 
thus expressed as 



+2ij(7f,a^^+af,7i^). 



kl 



(2.20) 



(2.21) 



where [ij] refers to the antisymmctrisation of indices such 
that for an antisymmetric tensor Xij = X[ij\ (and simi- 
larly later on, (ij) refers to the symmctrisation of indices 



such that for a symmetric tensor = X(y)). These 
Weyl components are explicitely given by 



C" 



Oi 



3t2 V^ + ^-S?')'^^' (2.22) 



(2.23) 
(2.24) 



and it is easy to notice explicitely that these two types of 
components are not independent and are related thanks 
to Eq. dMl)- We thus obtain that = Cf^^p^C'^'P'' is 
given by 



12C^ sin^ 



' 3VJ TT 

T+4 

= 3'CVi{l~Pi) 

= 3C^l + 3qif{2-3qi) 

= 6C^l + 2Qif {l-Qi) (2.25) 

where in the last equalities, / can take any of the values 
1, 2, 3. The function is given by [using Eq. ([^1^ ] 



1 



(2.26) 



The Weyl tensor and his algebraic properties are used 
to classify spacetimes by their Petrov type. Defining the 
tracefree symmetric rank-2 tensor Qij = —Eij — iHij 
in terms of the electric {Eij) and magnetic (-ffy ) parts 
of the Weyl tensor, provides a way to classify the Weyl 
tensors [1^ . For a Kasner spacetime, since the magnetic 
part of the Weyl tensor vanishes we obtain 



Oi 



Oj 



P,(1~P,)S}. (2.27) 



The eigenvalues A; of Qj define its Petrov type. They are 
all real and different so that a generic Kasner spacetime 
(tu 7^ f 7 f ) h^-s Petrov type I while the two axially sym- 
metric solutions with degenerate eigenvalues are of type 
O when = | and of type D when to = ^ [l^ . 



III. PERTURBATION THEORY IN A KASNER 
SPACETIME 

The Kasner spacetime being a particular case of 
Bianchi universe, the study of the evolution of the pertur- 
bations derives simply from our previous formalism P, Q 
that we specialize to the vacuum, using that S" /S = 
-1/(477^) so that n' = -2%"^ and the property ((^1^ . 



A. Generalities 

Following the formalism developed in Refs. [H, 0, we 
consider the general metric of an almost Bianchi I space- 
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time, 

ds^ = S^[-{l + 2A)df]'^ + 2Bidx'dT] 

+ (7„- +/iy)dxMx^] . (3.1) 

Bi and hij can be further decomposed into scalar, vector 
and tensor components as 

^ d^B + Br, (3.2) 
/ly = 2C(7.y +Sy) + 2a»9ji; + 2%£;j)+2^,j-, 

with 

5,5' = = d,E\ El = = d,E'^. (3.3) 

We then construct the gauge invariant quantities and de- 
fine the conformal Newtonian gauge by the conditions 

B = B' =E = 0, (3.4) 

so that 

A = <^, C = -*, <P' = -{E'y. (3.5) 
We also introduce the extremely useful variable [2^ 



X = ^ + ^ 



n 



(3.6) 



The tensor variable Eij is readily gauge invariant. At 
this stage, we are left with 10 — 4 = 6 degrees of freedom: 
two scalars ($ and \1/), two vectors ($*) and two ten- 
sors (Eij). Contrary to the perturbation theory around 
a spatially homogeneous and isotropic spacetime [l^l , the 
three types of perturbations do not decouple. 

As we shall now see, the perturbation equations will 
exhibit four constraints so that only two dynamical de- 
grees of freedom related to the gravity waves remain, as 
in Minkowski or de Sitter spacetimes. 



B. Mode decomposition 

The perturbation equations are conveniently written in 
Fourier space and we decompose any quantity in Fourier 
modes as follows. Using the Cartesian comoving coor- 
dinates system {x'} on the constant time hypersurfaces, 
we decompose any scalar function as 



The comoving wave co- vectors ki are constant, k[ ~ 0. 
We now define /c* = 7'-'fcj which is now time-dependent 
and cxplicitely given by 



h ( ^ 



-2Q, 



(3.8) 



so that 



(3.9) 



For any mode fc^, a basis {ei, 62} of the subspace perpen- 
dicular to ki can be constructed from the natural Carte- 
sian basis as (sec Appendix A of Ref. Q for the details 
of the construction) 



with 



cos 7 cos /3 cos a — sin a sin 7 
cos 7 sin a + cos a cos (3 sin 7 
— cos a sin /3 



— cos 7 cos (3 sin a — cos a sin 7 
cos 7 cos a — sin a cos (3 sin 7 
sin a sin (3 



(3.10) 



(3.11) 



(3.12) 



The three Euler angles (a, /3, 7) depend on time and are 
explicitcly given by 

sin7sin/3=— , cos7sin/3 = — (3.13) 
k k 

so that tan7 = X2/X1 = {k2/ki){r]/r]oY'^^^'^^'> and 

cos/3 = ^. (3.14) 
k 

The previous relations determine f3{r]) and 7(77). The re- 
quirement that {ei, 62} remains orthogonal to ki imposes 
that 



a' = — cos Pj' , 
which can be rewritten as 

, (Q2 - Qi) X3 1 



k + 

X2 X\ 



(3.15) 



(3.16) 



which determines a (77) up to an integration constant. 

The vector and tensor modes can then be decomposed 
respectively as 



$,(fc„77) = ^aih,v)enh), (3.17) 



(3.7) and 



(3.18) 



A= + ,x 

where the polarisation tensors have been defined as 



pipi _ p2 2 12 I 21 



V2 



V2 



(3.19) 
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C. Shear components 

The perturbation equations in Fourier space involve 
the components of the decomposition of the shear on the 
basis {k, ei, 62}. Since it is a symmetric trace-free tensor, 
it can be decomposed as 



A=+,x 



(3.20) 



This decomposition involves 5 independent components 
of the shear in a basis adapted to the wavenumber hi. 
We must stress however that {a^^ ^a^a^o-^x) must not be 
interpreted as the Fourier components of the shear, even 
if they cxplicitely depend on ki . This dependence arises 
from the local anisotropy of space. A similar decomposi- 
ton for Tiij defines the coefficients (E|| , S^ai ^t\) ^^^d we 
obtain from Einstein equation in vacuum that it must 
satisfy the constraint 



6 = = 



2 _ f^ij^ 



= |l]^+2^E^_„+^S^,, (3.21) 



which is independent of ki. 

The equation of evolution of the shear (|2.17p implies 
that 



11 va 



(3.22) 



b.X 



(3.24) 



a.b 

where the matrix M^b defined by 

The equation p.20p can be inverted to get 



2Q, 



(3.25) 

(3.26) 
(3.27) 

(3.28) 



The last equation implies, using Eq. p.l9|) . that 

i 

i 

D. Perturbation equations 



(3.29) 
(3.30) 



We use the equations derived in Refs. (T], Q when ap- 
plied to the particular case of a vacuum solution (so that 
V -> 0). 

The two scalar perturbations are cxplicitely obtained 
from the tensor modes as 



We then deduce that ^' is given by 



(3.31) 



-nx' 

2-E„ 



(3.32) 



' a,b,X 



while the second Bardeen potential is then given by 



A:2$ = /c^ (4 - E J ^' -f 3HX'. 



The vector mode is then given by 



(3.33) 



b II 



(3.34) 

This shows that the scalar and vector modes arc ob- 
tained algebraically from the tensor modes, which are 
the only degrees of freedom that propagate. Using the 
shorthand notation (1 — A) for the opposite polarisation 
of A, i.e. meaning that if A = +, then (1 — A) = x, and 
vice- versa, and introducing 



Ma = SE\, 



(3.35) 



we have 
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2-E„ 



^l,+HT.'ax. (3.36) 



The general solution of this linear equation can always 
be obtained in terms of four transfer functions as 

Ex{h,ri) = Tly{ki,i])ay{ki, Tjo) + Tj^y{ki,r])by{h,rio), 

(3.37) 

where ax' {ki, rjo) and b\i (ki, rjo) are the initial conditions. 
The notations J, and f refer respectively to a decaying 
and a growing mode. When the two polarisations are 
decoupled (which is the case asymptotically at early and 
late times, as we whall see below) then we only have two 
transfer functions to consider T^^ = T^Swi and T^-,^, = 
T^Sw since the transfer of power from one polarisation 
to the other is negligible. Similarly we define transfer 
functions for ^P, X and by 

'^(h,v) = Tl,^{ki,r])ax{ki,'qo) + T^.^(k„rj)bx(h,r]o), 

(3.38) 

X{ki,T]) = T]^^^{ki,r])ax(h,T]o) + T^^^{ki,r])bx{h,r]o), 

(3.39) 

'^a{ki,r]) = Tl^-^{h,T])ax{h,rjo) + T^^^^{h,ri)bx{ki,7]Q). 

(3.40) 



where is the background value of given in 
Eq. (|2.25p such that at the background level S = 1, and 
at the perturbed level up to second order we have when 

JWC^ = C^sW, (4.2) 
5(2)^2 = (4.3) 

The perturbation of S at first order in the perturbations, 
S'^^)(x, 77), can be decomposed in Fourier modes as 

SW(X,^)= /^sW(k,^)e'>^- (4.4) 

In order to assess the behaviour of at early and late 
time, we assume that at an arbitrary initial time rjo the 
initial conditions for the two modes defined in Eq. p.37p 
are such that (1) they are not correlated, 

{axik)b\,ik'))=0, (4.5) 

and that (2) they have the same initial power spectrum, 

(aA(k)a^,(k')) = <53(k - k')6xx'Pinit{h), (4.6) 



E. Summary 

In conclusion, for a given model, one can determine 
(S|| , T,^a, S^a) for each mode ki and then solve Eq. p.36p 
for the gravitational waves. One can then deduce $, ^ 
and algebraically. As expected, only two degrees of 
freedom can propagate. The existence of the vector and 
tensor modes arises from the fact that isotropy is violated 
so that the SVT modes do not decouple. The Kasner 
case, being a vacuum solution, is thus simpler than the 
Bianchi I case studied in Refs. [H, [I]- 

IV. STABILITY ANALYSIS 

A. Generalities 

In full generality, the stability analysis requires first to 
solve Eq. p.36p to determine the four transfer functions 
defined in Eq. p.37p . This allows to evaluate the effect 
of the perturbations on the square of the Weyl tensor, 
which is by construction independent of the choices of 
gauge and coordinates system. We thus introduce 



{bxik)b\,ik')) = 6'{k-k')dxX'Pnut{h). (4.7) 

The initial power spectrum is an unknown function of 
the comoving wave- vector ki. It is clear that 

(sW(x,77))-0 (4.8) 

at all time. This is indeed not the case for the pertur- 
bation of S at second order in perturbations, S^^^(x, 77), 
since it is quadratic. The previous definitions allow to 
compute that at lowest order (S(x, 77)) = (S(^)(x, 77)) 
with 

(S(2)(x,77)) = I |!^Pi„i,(fe0^sf (k,^) (4.9) 
and we define 

S(2)(k,77)^^4^)(k,7y). (4.10) 

A 

The function E\ (k, 77) can in principle be expressed in 
terms of the transfer functions that appear in Eq. p.37p 
and in terms of the coefficients E||(k, 77)... 

From a numerical point of view, the situation is thus 
clear but very time consuming. It is also unnecessary 
to solve the evolution equations in their full generality 
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since we are interested in the asymptotic behaviour of 
(S^^) (x, 77)) when rij-qQ — ^ +00 or 77/770 0. In these two 
regimes, the wave numbers tend to focus along a principal 
axis since [see Eq. 
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k ^ ki 



and 



k fca 



'70 



-Q3 



0. 



(4.11) 



(4.12) 



We shall thus discuss the dynamics of the perturbations 
when their mode is aligned along a principal axis (? IIVB[) 
and their contribution to S^'^-' (k, 77) IIV Cp . These two 
steps can be performed completely analytically so that 
we can then discuss the general behaviour of (S(^)(x, 7/)) 
at late time (? |IVD| and early time f? llVE)) . 



B. Dynamics of the perturbation for modes aligned 
along a principal axis 

1. General behaviour 

We assume that only one of the components of the 
Fourier space, indexed by /, satisfies kj ^ 0. It is then 
clear that fc^ = kj{r]/rio)^^'-^' so that 



2Qi 



(4.13) 



Then, since Sj^ =0, we conclude from Eq. p.22p that 
Y,a = because of Eq. ([5:^ and thus 



Sva = 0. 

This implies, from Eq. that 
$a =0. 



(4.14) 



(4.15) 



It follows from Eq. p.24p that E^^^ = 0, and it can be 
checked that 







where we have defined 

A/ = qma.K{i^I) — 9min(i^/) 

which is explicitely given by 



(4.16) 



(4.17) 



A, = A |cosn7,| = y^-3g| = 2^^1^.(4.18) 
The scalar perturbations are then given algebraically by 



X 



2(1 -Q/) 



(4.19) 



3 -i-:: 



FIG. 2: Evolution of the index of the Bessel function, , 
with vj (v^ in solid line, i^^ in dashed line and in dotted 
line). We recall that vf — whatever m. 



that derives from Eq. p.3ip from which we deduce that 
Eqs. (|3.32II3.33P reduce to 

k^^ ^ -nx', k^<^ = -{l-2Qi)HX' .{A.2Q) 



The only equations to solve are Eqs. p.36p for the grav- 
ity. They decouple for each polarisation and lead to the 
system 



-■2Q1 



1 

4^2 



ki[^ 

Vo 



-2(3/ 



+ 



1_9A| 
4772 



0, 



(4.21) 



=0.(4.22) 



These two equations are compatible with the general 



equations 



obtained for w = ^ and, in this 



particular case, with those derived in Ref. [J]. The gen- 
eral solutions of Eqs. (j4. 21114.22]) can both be obtained as 
[see Eqs. (jCHIC2p ] the linear combination (|3.37p with 

T2(fc/, 77) = J,] {ki,r^), Ti{ki,r^) = N,) {ki, 77) (4.23) 
where 

kirjo 



Ji\ki,ri) = J, 



1-0/ 
1-0/ 



V_ 

Vo 



1-Qi 



(4.24) 



(4.25) 



Here, J^a and iV^A are the Bessel and Newmann functions 



of index 



3A7 



2(1-0/) 



/3(l + 0/) 

(1-0/) ' 



(4.26) 



and is depicted on Fig. [5] (note that i/j is always non- 
negative while can be negative). It is then clear that 



\X\ 



V2 



(4.27) 
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so that the two Bardecn potentials are given by 

-2(1-Q/ 



* = -- 



1 



V 



$ = (l-2Q/)«' 



^T+ dE+ 
4(1 -Q,)mniL ^4.28) 

(4.29) 



If follows that for these particular modes, the solution 
of the dynamics is obtained completely analytically. Fig- 
ures [3] give an example of the behaviour of the tensor and 
scalar modes. 



2. Late time behaviour 



with 



1 



(4.33) 



These expressions hold as long as ^ 1, that is for 
tn 7^ ^, a particular ease that is discussed in Appendix IB] 
The transfer functions for the + polarisation behave as 



Tl 



Ti 




2(1-Q3) 



1-Q3 



(4.34) 
(4.35) 



The late time behaviour of the perturbations can 
be easily understood using the large argument expan- 
sions (IC3I) and (IC5I) of the Bessel functions, since 



with 



nil \r]o 
V_ 

nil \t]o 



i-Qi 



sm 



'70 



i-Qi 



2kir]o ki7]o 



2-3(71 



1 - I 



(4.30) 



It follows that the late time behaviour is in general of the 
form 



V TT^l \7]o 



1-Ql 



a[°°^ cos 



1-Qi 



(00) 
(4.31) 

where and arc two constants that can easily be 
obtained in terms of a\ and bx- X is then proportional 
to E+ [see Eq. (|4.27p ] and $a = 0, which completely 
determines the solution [see Eqs. ()4.28M.29|) ]. 



3. Early time behaviour 

Asymptotically, the dominant term of the wavevector 
is aligned along the direction expanding the fastest, i.e 
1 = 3. The two polarisations behave differently because 
A3 7^ 0. Using the small argument expansions (|C4[) 
and (|C7p of the Bessel function, we obtain that, since 
1 - Q3 > 0, when 7] ^ 



V 



1 



r(i + z.3> 
r(^3") 



2 

m 



IL 

1-Q3 



(4.32) 



C. Weyl tensor evolution for modes aligned along a 
principal axis 

For these modes, the expressions of 3^^^ (k, 77) and 

Sj^ (k, 77) turn out to simplify greatly and can be ex- 
pressed in terms of the dimensionless parameter 

X = kr] = xir/ . (4.36) 

First, S^^'' (k, rj) can be expressed in terms of the two 
transfer functions and their first derivatives as 



1 V fW^!Zi_ 

2(1 - Qi)il + 2Qi)x^ ^^/+ dilnr^y 



The two coefficients f^^ can be computed and are given 
by 



(0) 

(1) 



V2Ai[-6{l + 2Qi)x^ + 2x% (4.37) 
V2Ai[-9{l + 2Qi) + (7 + 2Qi)x^]. (4.38) 



Interestingly, the x polarisation does not contribute. In 
the above expression T+ stands either for t| or t| de- 
pending on the solution (decaying, growing, or any linear 
combination of them) that is considered. When A/ = 0, 
"Ey^^ vanishes identically and, again, since (S*^^)) = we 
shall not consider this linear part any further. 
The second order S^^^ can be expanded as 



1 



6(1 + 2Q/)2(1-Q,) 



(4.39) 



EE/. 

1=0,1 j=0,l 



9(ln77)* d(\x\rj)3 



with the same conventions as for The eight coeffi- 

(i i) 

cients are explicitely given by 

= 3(1 + 2Q/)2(13 + 8Q/)(1 + Q/) 

+6x^{-l + 4Qj) +8x^ (4.40) 
/(^o-") ^ 6(1 + 2Q,)2(l-g2) + 18^2 -K8a;^(4.41) 



9 




f(0,l) 



f(0,l) 
J X 

and 

(1,1) 



^(1,0) 

\ [3(1 + 2Qif{l + Q/)(22 + 25Q/ - 2Q]) 

+2(1 + 2Q/)(-26 - 29Q/ + Q])x^ 

+4(1 - Q,)a;4] (4.42) 

(1,0) 



= /: 

= -A{l-Qi)[{l + 2Qif ^x"] 



(4.43) 



/- 



[9(1 + 2Qj)2(l + Qj)(3 + 2Q/ + AQ]) 



J X 



(1,1) 



X (5 - 2Q/) + 2(13 + 4Q/ - 8Q?)a;'' + Sx^ 
+9(1 + 2Q7)(1 + Q/)(-17 - 18Q/ + SQ^x^] 

= 2[{l + 2Qif +Ax^]. {AAA) 
D. Asymptotic behaviour at late times 



At late time {"q — > oo), only the growing mode domi- 
nates and we can safely neglect the decaying mode since 
it would only lead to a redefinition of the phase (/9oo ; see 
Eq. (|4.3ip . As can be concluded from the evolution (j4.1ip 
of the wave number, the parameter x defined in Eq. (|4.36p 
behaves as 



i-Qi 



(4.45) 



which is always a growing function of rj (remember gi < 
and thus Qi < 0; see Fig. [1]). We thus need to take the 
limit when a; oo of the transfer functions and of 
their derivatives. From Eq. (j4.3ip . we conclude that 

(4.46) 



,/2(i-gi) 



1 - 



Since 



dx 
3 In 77 



{1-Qi)x, 



(4.47) 



it implies that, at leading order in x. 



dhn] 



2xil-Qi) 



sin 



1 -Oi 



(4.48) 



Using the expansion (|4.39p . we conclude that 



Ax^ 


\ain7jj 




3 


(1 + 2Q 


i)2(l-0i) 



37r(l + 2Qi)2 



(4.49) 



This expression is valid for any Kasner spacetime but 
for the particular case nj = ^ that is discussed in Ap- 
pendix |B] (since in that case Qi = — ^). 
It follows that 

E^P(k,T])(xx^ (XT]^^^-'^'\ (4.50) 

which is unbounded at late time. 

E. Asymptotic behaviour at early times 



As can be concluded from the evolution (|4.12p of the 
wave number, the parameter x defined in Eq. (|4.36p be- 
haves, when 7] — > 0, as 



1-Q3 



(4.51) 



which is always a decreasing function of rj (except if = 
^ , see Appendix [B] for that particular case) . It follows 
that we need to evaluate the limit a; — >■ of the transfer 
functions. 

At early time, one can however not disregard the de- 
caying mode since it diverges faster when ry — >■ 0. The 
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expression (|4.39|) behaves as 



Discussion 



(4.52) 



(1 - Q3)x^ 

3(1 + Q3) 
2il-Q3)xi 



(22 + 25Q3 - 2Qi)T+ 
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(5-2Q,)(3 + 2Q3+4Q^ 



4 9T. 



1 



9 In 77 
9Tx 



3 "ainr; 3(1 -Q3) V^lnr; 



In order to determine the behaviour of the transfer func- 
tions, we consider two cases. 

First, if we consider only the growing mode [i.e. a\ 7^ 
and = in Eq. p.37p ] then we deduce from the 
expressions (|4.32p and (|4.35p of the transfer functions 
that 



-A3TI 



dTl 

dTl a;^ 

9h^ " "2(1 -Q3) 



(4.53) 



At leading order, the contributions of the two polarisa- 
tions to S^^) are thus 



During the cosmic evolution, for any modes, the com- 
ponents of the shear evolve from an early value corre- 
sponding to the situation of the mode being aligned to 
the third spatial axis (index / = 3) to a situation corre- 
sponding to the mode being aligned to the first spatial 
axis (index / = 1). This dynamics arises from the fact 
that the Euler angles are not constant over time and that 
/3 evolves from to ^ and 7 from ^ to (except in the 
axisymmetric cases ■07=^,^ for which the value of 7 
does not matter). This implies that Yi^^ varies from iq^ 
to 3^1, and S^x from to but with a transient de- 
viation from and from -^^3, to "^^ii as depicted 
on Fig. m 

It follows that the asymptotic behaviour of the pertur- 
bation of the Weyl tensor can be discussed by looking at 
these two particular regimes. Figure [5] shows how a nu- 
merical solution interpolates between these two regimes. 

Using this procedure, analog to the one developed in 
Ref. U , we conclude that the Kasner spacetime is unsta- 
ble at late time since J'^^^C^ decreases with time slower 
than C*^, which means that the Weyl tensor of the per- 
turbation always dominate the Weyl tensor of the homo- 
geneous space at late time. Indeed, for a given mode the 
late time asymptotic regime is reached when 



9(l-4g3)'(l + Q3)yt2 



8(1 -gs 



4')(fc,,77) ^ 2(1 + Q3- A3)T;f 



(4.54) 



Eqs. and (|i35|) tells that - x^^ and - x° 

so that S^^\ki,ri) -> (since > 0) and E''^\ki,rj) 
const. Thus, the contribution of the growing mode leads 
to a bounded contribution to S'^^. 

On the other hand, if we consider the decaying mode 
[i.e. ax = and bx ^ in Eq. p.37p ]. i.e. the most 
diverging mode when x — >■ 0, then 



dlni] 

dTj ^ 2(1 - Q3) 
9 In 77 TT 



(4.55) 



We conclude that at leading order the contributions of 
the two polarisations to S^^^ are 

(k,r7) ^ ^il_Ml(5- 203)(3 + 2Q3 + 4Ql) 



S(f^(k,77) -> 2(1 + Q3 + A3)ri 



3; 
(4.56) 



Since, from Eqs. and (|i35|) . (x x'''^ with > 

3 for all cases, the x contribution diverges faster when 
77 — >■ and we conclude that 



(4.57) 



77 > 77jf*'' = max 



k2\ (Q2-Q1) 



ki 



(Q3-Q1) 



(4.58) 

At any given time 77, we denote by S(k, 77) the set of 
modes such that r] > r]^*^'^ so that Eq. (|4.9p implies that 



d^k 



(S^^^(x,77)) > / _^Pi„i,(fc,)s(2)(k,r;), (4.59) 



i2n) 



since S^^^ (k, rj) is positive. Given that S^^^ (k, 77) is 
converging toward 5(^)(fci,77) which is diverging at late 
times, we deduce that S^^^ (k, 7;) is diverging at late times. 
This is compatible with the analysis of Ref. Q for the 
case of an axisymmetric Kasner spacetime, as well as 
with earlier findings 0, Q on the amplification of gravity 
waves in anisotropic inflation during the shear dominated 
era, which can be described by a Kasner era of finite du- 
ration. 

At early time the conclusion depends on the choice of 
the modes that are considered. We have shown that if 
decaying modes are present then the spacetime is unsta- 
ble when 77 — >■ 0, following the same argument as above 
but with the modes such that 



^, early 
77 < 7/^ 



(Q3-Q2) 



\ (Q3-Q1) 

(4.60) 

whereas it is stable if one imposes to have only growing 
modes since then S^^-* remains bounded. As pointed out 
in Q (see Appendix B), this result is compatible with 
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the analysis by Ref. [g where it was pointed out that 
a condition on the perturbations should be imposed for 
the spacetime to be stable when 77 — > 0. As explained in 
Ref. [4j , this condition precisely kills the decaying modes 
A/"/ for both perturbations. 

V. AXIALLY SYMMETRIC CASE WITH = f 

In order to compare our analysis to the existing liter- 
ature, we consider the particular case of an axially sym- 
metric case with w — ^. This particular case was studied 
in Ref. Q and corresponds to the axisymmetric case 

2 1 , ^ 

91 = -g, 92 =173=2 (5.1) 

so that 




-2-10 1 2 



Log(;7/;7o) 

FIG. 4: Evolution of E II (solid), E^i (dashed) and (dot- 
ted) for k\r]o = fcxi^o = 1 (heavy) and fciryo = 10, fcxr;o = 1 
(light) or k-iTjo — 1, k±i]o — 10 (light) for a Kasner universe 
with = 



Ai = 0, 



A, = A, = 1. 



(5.2) 



In that case, and in that case only, we perform a per- 
mutation of the axis (3 — >■ 1, 2 — >■ 3 and 1 — >■ 2) with 
respect to the parameterisation of angles given in § IIII Bl 
so that the special direction is the first axis (labelled by 
the index / = 1). 

After deriving the gravity wave equation (§ IV A[) and 
showing that the two polarisations are decoupled at 
all times (and not only asymptotically as for a generic 
case), and expressing the vector and scalar perturbation 
(§|TE|), we investigate the late (§|VC]) and early (§|VD| 
time behaviours (and the interpolation of the exact nu- 
merical integration between these two asympotic regimes 
in §|VEl). 



A. Gravity wave equation 



We set k'j^ = fc| + fc| , so that Eq. (|3.9p takes the form 

e = kj(^)\ki(^y' (5.3) 

and we shall assume that kik± ^ 0, since otherwise we 
are back to § IIVBI Because of the rotational invariance, 
Eq. p.l4|) implies that the Euler angle 7 is constant 



tan 7 



(5.4) 



so that Eq. (|3.15p implies that a' = 0. We can thus 
choose 

a = 0. (5.5) 
The Euler angle /? is obtained from Eq. (|3.13p 



cos^=— — , sm/?=— — 
k \r]o J k Vr/o 



(5.6) 



It follows that Eqs. p.im3.12p simplify to 

cos 7 cos /? 



Wi" = I COS /3 sin 7 

— sin [3 




(5.7) 



The components of the shear then take very simple ex- 
pressions. Eq. p.26p reduces to 



while Eqs. 



sin^ ;3 - 2 cos^ ^, 



give 



1 = 3 sin /? cos /?, 
and Eqs. p.29ll3.30|) lead to 
3 2 



0. 



(5.8) 



(5.9) 



(5.10) 



These expressions satisfy indeed Eq. (|3.2ip and are de- 
picted on Fig. m We recover the behaviour (|4.13p for , 
(liTi]) for and for S^a- 

Since I].px = 0, the two polarisations decouple and 
their equations of evolution p.36p simplify to 



Ai'x.-H + Wx,-H('=^>Ai)/^x,+ = 



with 



2 

T + 



S^l)] 



1 



1 + 21]^ ' 



v ^ 



(5.11) 



(5.12) 



(5.13) 



where we have used that E'^ = — 2HI]^]^ and Ef^^ = 

\f2'KY?^^. We emphasize that our expressions agree with 
those of Ref. Q (see for instance Eq. (40) of that refer- 
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B. Scalar and vector perturbations 



and 



The expressions of the scalar perturbations are there- 
fore 



X 



2-S„ 



(5.14) 



that derives from Eq. p.3ip from which we deduce that 

3 




(5.15) 



-B4 



and 



V2- 



(5.16) 



The absolute value of the coefficient in Eq. (|5.14p evolves 
from /V2 to i^f /\/2 = 3/\/2 and $a vanishes 
both at early and late times. 



C. Late time behaviour 

At late times, the wave-number behaves as fc ~ kir]/r]o 
so that X ^ fci??o('y/%)^- For the two polarisations, the 
growing and decaying modes entering the solution p.37p 
behave as 



Jn 



kiVo f V_ 
2 V'?o 



Jo (I) (5.17) 



No 



(5.18) 



The two modes combine to modify the phase so that the 
general solution is 



Ex 



2A 



(oo) 



(oo) 



(5.19) 



The behaviour of the scalar and vector perturbations are 
therefore easily obtained since /3 — )- f , i.e. 



sin /3 



cos/3 ^ 1, 



(5.20) 



so that ~ —2 while E^i ~ 3sin/3 and S^- 
-3sin2^/V2. It follows that 



X 



-3 (''A (v 



E^ 



(5.21) 



-3i^/2 /fc_L 
hm V ki 



-7/2 



E^ 
E^ 



(5.22) 



These expressions can then be inserted in Eq. (|4.39p that 

is valid for tu = ? to conclude that at late time, the 
6 ' 

polarisations contribute in the same way to S'^^ as 



:(2) 
"A 



8x^ 
"37 



(5.23) 



which is equivalent to our prior result (|4.49|) . The result 
for the X polarisation matches the results of Ref. [1] up 
to normalisation factors of the transfer functions, but it 
appears to be different from the conclusions on the + po- 
larisation. Asymptotically, the Fourier mode is aligned 
with the preferred direction of the spacetime, and it turns 
out to be difficult to compare our result with Ref. [3| 
given that most expressions of their formalism are sin- 
gular for a mode along that special direction. We are 
confident that our result is correct in that limit since 
from Eq. (|4.49p we find that C^S^^^ which is the per- 
turbation of the square of the Weyl, when expressed in 
terms of the transfer function takes the form 



2^(2) 



— > 



54 



dTl 



k^T 



(5.24) 



This is formally the same expression as for a pure 
Minkowski space-time [see (Eq. IB5[) ]. which is expected 
when one takes the small scale limit. 



D. Early time behaviour 

At early times, the wave-number behaves as ~ 
fcj^ (77/7/0 )~ 2 so that X ~ k±r]o{r]/r]o)2 ^ which goes to zero. 
As in § IIVEI we cannot neglect the decaying mode. We 
can use the expression of S(2) obtained in § HVEl since it 
is well-behaved when 93 — > ^ . 

If we first consider the growing mode [i.e. a\ ^ and 
6a = in Eq. p.37p ] then we deduce that 



Tl 



Jn 



J3 



2A;±77o — 



2k±T]o 



from which it follows that 
dTl 



3„t 
2 



dTl 



M2x) 
Jsi^x), 



2rr\ 

-X 



dlnr] 2 ^ dim] 
Asymptotically, wc have from Eq. ()C4p that 

„3 



Tl 



1 



T 



X 

~6 



(5.25) 
(5.26) 

(5.27) 
(5.28) 
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Since 



VI. CONCLUSION 



cos /? ' 



kj_ \rio 



(5.29) 



S|| ~ 1 while ~ 3cos/3 and ~ —3/^/2 and it 
follows that Tx,x = 7*0, x = and 



-3T4 



(5.30) 

With these expressions, we conclude that at leading or- 
der, the contributions of the two polarisations to S*^^) are 



-3iV2 fci f ii\ f 



27 



7(2) 



X 

36' 



(5.31) 



which both remain bounded. 

On the other hand, if we consider the decaying mode 
[i.e. flA = and h\ ^ in Eq. p.37p ]. i.e. the most 
diverging mode when a; — )■ 0, then 



m 



2fc_L?7o 



= iVo(2a;) (5.32) 
= iV3(2a;), (5.33) 



the asymptotic behaviours of which are 



TX In (x) 

n 



X 

n 



-3 



(5.34) 



It allows to conclude that the contributions of the two 
polarisations to E^^^ are 



90 



"(2) 



20 



(5.35) 



As expected both terms diverge when 77 — > 0. 



E. Interpolation between the two regimes 

The transition between the two asympotic regimes is 
illustrated on Fig. |4] for the components of the shear. 
While the transition is quite sharp, it is shifted depending 
on the value of k±/ki. 

Then, we can integrate numerically Eq. (|5.1ip for dif- 
ferent modes. This exact numerical solution is compared 
in Fig. [5] to the two asymptotic forms, which demon- 
strates that the transition is sharp even for the evolution 

of the modes. In Fig. [6l the full numerical solution of 

(2) 

Sj^ ' (k, 77) is compared to the two asymptotic behaviours 
described in S|VC] and WDI . This shows explicitely the 
validity of the asymptotic expansions. 



We have performed the stability study of a generic Kas- 
ner spacetime with respect to linearperturbations using 
the formalism developed in Refs. [if, Q. Since Kasner 
spacetimcs are solutions of the vacuum Einstein equa- 
tions, one expects only two degrees of freedom to prop- 
agate, which we explicitely demonstrate. They can be 
identified with the two polarisations of gravity waves. 
The extra scalar and vector components of the metric are 
then obtained algebraically from constraint equations. 
Contrary to a Minkowski or de Sitter spacetime, they 
do not strictly vanish because the spatial sections are 
not isotropic anymore. 

We have shown that for any Fourier mode, the be- 
haviour of the gravity waves interpolates between two 
asymptotic regimes, at early and late times, in which the 
mode can be considered almost aligned with a princi- 
pal axis of the spatial metric. Interestingly, in this limit 
the gravity wave equation can be integrated analytically. 
This allows to set a lower bound on the square of the Weyl 
tensor generated by the perturbations and compare it to 
the one of the background. We have concluded that at 
late time the Kasner spacetimes were unstable with re- 
spect to linear perturbations but in the particular case 
of n7 = ^ , which corresponds to an axisymmetric config- 
uration, product of a two-dimensional Milne spacetime 
and a two-dimensional Euclidean spacetime, which ac- 
tually maps to one quarter of the Minkowski spacetime. 
At early time, the conclusion depends on the modes in- 
cluded in the analysis. If one includes decaying modes, 
the Kasner spacetimes are again unstable while stability 
to linear perturbations is recovered only if the growing 
modes are excited. This latter result was already known 
since Ref. 0] and our analysis confirms the one of Ref. Q 
that was limited to the particular case of an axisymmetric 
configuration with nj = ^. This result is also compat- 
ible with the amplification of gravity waves during the 
shear dominated era of an anisotropic infiationary phase 
described by a Bianchi I spacetime [3, 0] . This is ampli- 
fication of a remnant of the instability described in this 
article for a finite Kasner era. 

This result has some importance concerning the dy- 
namics of the universe close to the singularity. On the 
one hand, the succession of Kasner era in the BKL for- 
malism and the rotation of the Kasner axes from one era 
to the other was studied by assuming homogeneous be- 
haviour, i.e. assuming the inhomogeneity scale was larger 
than the mean Hubble patch . The instability at early 
time (i.e. going toward the singularity) can alter the con- 
clusions on the approach of the singularity, even though 
each Kasner era has a finite duration. In particular, as 
soon as the effect of the gravity wave is large enough, the 
backreaction would have to be taken into account, and 
it is not clear that the evolution toward the singularity 
follows a series of BKL oscillations. This requires fur- 
ther investigations that go beyond perturbation theory, 
see e.g. Ref. [2^. On the other hand, the instability at 
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10-^ 0.001 0.01 0.1 1 10 100 lo-i 0.001 0.01 0.1 1 10 100 10-1 0.001 0.01 0.1 1 10 100 



V/lo V/Vo l/lo 

FIG. 5: Evolution of a mode with kirjo = k^rjo = 10 and fc2 = in a Kasner universe with vj = ^: (left), T| (middle) and 
TjJ (right) compared to their asymptotic late time (dashed) and early time (dotted) asymptotic forms. The agreement at early 
times is very good so that the continuous and dotted curves are hardly distinguishable. 




FIG. 6: Evolution of (k, ry) for a mode with fci?7o = ^3770 = 10 and ^2 = in a Kasner universe with 137=^ for the 
polarisation x (left) and + (right) as well as their late time (dashed) and early time (dotted) asymptotic forms. 



late time could also be relevant in models where the pre- 
inflationary era is described by an anisotropic era. The 
gravity waves modes generated during this era can po- 
tentially be observed [1, 0, [13] if the number of e-folds 
during inflation is not too large and their presence can 
modify the onset of inflation. 

Besides the speculative phenomenology of the early 
universe, our result sheds some light on the peculiarity of 
the Kasner spacetimes among the class of homogeneous 
vacuum solutions of Einstein equations. 
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Appendix A: Different parameterisations of Kasner 
exponants 

There exist several ways to parameterize, the coeffi- 
cients {pi,P2,P3) of the Kasner metric. Because of the 
two constraints (|2.4p . it is in fact a 1-parameter family of 
triplets that is completely specified by only one number. 
Three parameterisations have been widely used. 

• Single exponant parameterisation. As we have seen 
in § III Al it is sufficient to give one of the Kasner 



15 



exponants to specify the two others. For instance, 
fixing the smallest one, pi , gives that P2 = I —pi — 
Ps so that pI + (1 — Pi — Ps)^ +pl = 1, from which 
we deduce that 



1 



V(l-Pi)(l + 3pi) 



1 r 



P3 = 



1 -pi + V(l-Pi)(l + 3pi) 



(Al) 



(A2) 



The two constraints (|2.4p allow to deduce that 
the Pi satisfy the following useful relations. First, 
{J2iPif = f2iP'i so that 

P1P2 + P2P3 + Pd.Pi = 0, (A3) 

also conveniently rewritten as 

Y.p-'=Q. (A4) 

i 

From i^iPi)'^ — 1, one deduces that 

i i 

This tells us that 

- 1) =P2P3 (A6) 

which implies in particular p\ + p\ + P1P2 ~ ipi + 
P2) = 0, with similar relations obtained after per- 
mutation of the indices 1— >-2— >-3— >-l. In terms 
of the coefficient qt, it leads to the relations 



292.3 = -qi ± Al , 



Al 



3g? 



(A7) 



In terms of the Qi = 3/2qi these relations take the 
simple form 



2Q2,3 = -gi±v3(i-Q?). 



(A8) 



This shows that all the properties of the Kasncr 
spacctime, but also the evolution of its perturba- 
tions when a mode is aligned with an eigendircc- 
tion of the shear, arc characterized by the choice 
of a single Kasner exponent corresponding to that 
particular direction. 

• Angular parameterisation. As we have seen in 
§ III A| the constraints p.4|) imply that the triplet 
can be specified by a choice of an angle tu S [f , f ] 
and its expression is given by Eq. (|2.5p and we have 
seen that 



A/ = -^ICOSTU/I 



(A9) 



• BKL parameterisation. Refs. 
eterisation 



Pi = 



a{u)' 



P2 



1 + u 
cf{u) ' 



P3 



use the param- 
u{\ + u) 



a{u) 

(AlO) 

with a{u) = 1 + M + u^. It can be obtained l2l[ by 
solving the two constraints (|2.ip once having set 



P3 
P2 



(All) 



and the condition pi < P2 P3 imposes that u G 
[l,-|-cxj[, with u = I corresponding to zu = ^ and 
U +00 to ZD — ^. 



Appendix B: Axially symmetric case with zu = ^ 

The particular case uj = ^ requires special attention. 
It corresponds to a spacetime with metric 



-dt^ + t^dz^ + Ax^ +Ay^, 



(Bl) 



which appears to be the product of a two-dimensional 
Milne space with a two-dimensional Euclidean space. It 
is the Taub representation [2^ of the Minkowski metric. 
With the change of coordinates defined by 

(t, X, y, z) ^ (T ~ t coshz, X = x,Y ~ y, Z = t sinhz), 

(B2) 

this metric is rewritten as 



dz^ + dA:2 + dy^ 



(B3) 



i.e. has a Minkowski metric but only covers the patch 
~ > 0, the upper cone {T > 0) corresponding 
to an expanding universe while the lower one {T < 0) 
describes a contracting spacetime. This explains why 
= [Eq. (12:251) ]. 
The choice of the time slicing plays an important role 
when studying the perturbations. In Minkowskian co- 
ordinates the spacetime is explicitely homogeneous and 
isotropic so that the tensor modes decouple from the 
other types of perturbations and evolve according to 



(a| -A 



{X,Y,Z}, 



hij ^ 0. 



(B4) 



Their amplitude is thus constant and the perturbation of 
the square of the Wcyl tensor 5C^ remains constant. Un- 
der its Kasner form, this no more obvious and one needs 
to check the consistency with the result in Minkowskian 
coordinates 



drj 



(B5) 



where the two modes have constant amplitude (so that 
the distinction between growing and decaying is no more 
relevant). 
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1. Description of the perturbation from a Kasner 
point of view 

a. Gravity waves equation 

We follow the general formalism of this article to study 
the perturbation in the particular case 



1 2 

91=92 = --, ^3 = g 



so that 



Ai = A2 = 1, A3 = 0. 



It implies that 



(B6) 
(B7) 

(B8) 



Obviously, this case can be deduced from the case 137=^ 
by 9i — > —Qi so that the Euler angles are again 



tan 7 = — , a = 0. 

K2 



(B9) 



It follows that the expressions for the shear compo- 
nents are also given by Eqs. (|5.8ll5.10l) but with all signs 
changed, ie. E|| — >■ — E|| etc. The evolution of the gravity 
waves is dictated by Eq. (|5.1ip with 

= - sin^ /3 + 2 cos^ P, (BIO) 
E^i = -3 sin (3 cos /3, S^a = (Bll) 

and 

^ -^ain^ f3, S^x = 0. (B12) 
v2 

The only difference concerns the Euler angle (3 which is 
now given by 

cos/3 = M^)~\ sin/3. %(^)\ (B13) 
k \rioJ k \rio/ 

b. Late time behaviour of the perturbations 

In the limit rj —>■ +00, k ^ k±{r]/riQ)i so that x ^ 
ki_rj{){j]/ri{))^ . This implies that 



cos (3 



k± \r]o 



sin /3 ^ 1 



(B14) 



^ —1, ^ — 3cos/3 and ^ 3/^2 so that we 
conclude 



The solution of the gravity wave equation is then ex- 
plicitcly given by 



Ex 



Af^' cos 



(00) 



(B16) 



for the two polarisations, and where the growing and de- 
caying modes are combined together, leading to a redef- 
inition of the phase. 

Since = 0, we cannot use the variable 5 anymore 
and thus work directly with S^^^C^. In the limit 77 — >■ cx), 
we obtain 



8x^ 
24x3 



d In rj 



X Tl 



cx rj 



(B17) 



It follows that (5'^'C| decays at late time so that the 
spacetime is stable against the perturbations. Note that 
this is the same behaviour as for the case w ^ but the 
main difference arises from the fact that (7^ = since in 
the general case decreases with time but slower 

than so that the effect of the perturbations dominate 
over the background quantity at late time. The stability 
of the Taub metric (jBip was studied in Ref. [IBl with the 
same conclusions as our analysis. 



Appendix C: Bessel functions 

The differential equation 



7-n 



1 .2 2\ 1 



has the general solution 

u{x) ~ \fxZi, \iix'*\ , 



= (CI) 



(C2) 



where is a linear combination of a Bessel function of 
the first kind (J^^) and of the second kind (Ny or New- 
mann function). 

We have that in a; ~ cx), 



while in a; 0, 



X — V— 7 

2 4 



r(i + ^.) 

We have that in x ^ 00, 

N,{x) ^ 
while in a; 0, 

N^(x) 
N,ix) 



1 /xy 

V2/ 



2 4 



In [-) + 

T{v) f2Y 



(C3) 
(C4) 

(C5) 

(C6) 
(C7) 
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